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SELF-SIMILAR SOLUTIONS TO THE CURVE 
SHORTENING FLOW 

HOESKULDUR P. HALLDORSSON 



^ ' Abstract. We give a classification of all self-similar solutions to 

the curve shortening flow in the plane. 
ON ' 

O ' 1. Introduction 

Q 

Some examples are known of immersed curves in the plane which 
move in a self-similar manner under the curve shortening flow (CSF). 
The two obvious ones are the straight line, which is not affected by 
the flow, and the circle, which shrinks homothetically to a point in 
finite time. The circle is the only simple closed curve which shrinks 
homothetically under the flow but there also exist immersed closed 
curves with the same property. They were classified by Abresch and 
r-~ Langer in |AL] . A discussion about homothetically expanding curves 

O ! can be found in [I] and [U]. The Grim Reaper curve is the graph of 

the function — logcos(x). Calabi discovered that it translates upwards 
with constant speed under the flow. In [A] , Altschuler shows a picture 
of a yin-yang spiral and claims that it rotates under the flow. 

So we have curves which shrink, expand, rotate and translate under 
the flow. But are there any other self-similar solutions, which combine 
two or more of these basic motions? The answer to that question is yes 
and in this paper we find and classify all immersed curves in the plane 
which move in a self-similar manner under the CSF. Our main results 
can be seen in the following list. 
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Translate: Only the Grim Reaper curve; Figure [U 

Expand: A one-dimensional family of curves. Each is embedded 
and asymptotic to the boundary of a cone; Figures [2J12J 

Shrink: A one- dimensional family of curves. Each is contained 
in an annulus and consists of identical excursions between its 
two boundaries; Figures HH3 

Rotate: A one-dimensional family. Each is embedded and spirals 
out to infinity; Figures IEHHJ 

Rotate and expand: A two-dimensional family of curves. Each 
is embedded and spirals out to infinity; Figures [TUirTTl 

Rotate and shrink: A two-dimensional family of curves. Each 
has one end asymptotic to a circle and the other is either as- 
ymptotic to the same circle or spirals out to infinity; Figures 



In Section [2J we show how the problem of finding all self-similar 
solutions to the flow can be reduced to the study of a two-dimensional 
nonlinear system of ODEs. The system has two parameters, A and B, 
which determine the type of the self-similar motion. The values of A 
and B are divided into four cases which yield curves with significantly 
different behaviour under the flow. In Sections [3H6] we look at each case 
separately and prove some properties of the corresponding curves by 
studying the solutions to the system of ODEs. 

Section [7] contains figures of the curves. 

2. All self-similar solutions found 

Let X : R — y C be an immersed curve in the plane. A self-similar 
motion of the curve is a complex valued function X on R x I of the 
form 

(1) X{u,t) = g{t)e if{t) X{u) + H(t) 

where I is an interval containing and /, g : / — y R and H : / — y C are 

different iable functions s.t. /(0) = 0, g(0) = 1 and H(0) = and hence 
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X(u, 0) = X(u). The function / determines the rotation, g determines 
the scaling and H is the translation term. 

By definition, this motion is the curve shortening flow of the curve 
X (up to tangential diffemorphisms) if and only if the equation 

dX 

( — (u,t),N(u,t)) = k(u,t) 

holds for all (u, t) G R X J. Here N(u, t) = i T(u, t) is the leftward 
pointing normal to the curve u i-> X(u, t) and the (signed) curvature k 
is given by k(u,t) = («(«,£), N(u,t)), where k is the curvature vector 
of the curve. Straightforward calculations yield that this equation is 
equivalent to 

g 2 (t)f'(t)(X(u),T(u)) + g(t)g'(t){X(u),N(u)) 

(2) +g(t)(e- if VH'(t),N(u)) = k(u). 

Now, we want this equation to hold for all (it, ()gRxI. In particular, 
it must hold when t = 0. Thus, the original curve X has to satisfy the 
equation 

(3) A(X(u),T(u)) + B(X{u), N(u)) + (C, N(u)) = k{u) 
for all «GR, where A = f (0), B = g'(0) and C = H'(0). 

First we look at the case where the translation term vanishes so 
C = 0. Assume X is an immersed curve satisfying the equation 

(4) A(X,T)+B(X,N)=k. 

Now, if we pick the functions / and g such that g 2 (t)f'(t) = A and 
g(t)g'(t) = B for all t & I, then it is clear that equation (J5J) holds for 
all (it, t) & R x I. A valid choice for / and g is 

,£log(2Bt + l) if 5^0, , 

/(*) = { 2B K ' ^ and g(t) = y/2Bt+l. 

At if B = 0, 
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Hence, under the CSF, the curve X moves in the self-similar manner 
governed by / and g as shown in equation (JT]). In other words, it ro- 
tates around the origin (unless A = 0) and dilates outwards if B > 
and inwards if B < 0. This motion will be called screw-dilation. 



What happens if we include the translation term HI Assume that 
the curve X satisfies equation ([3]) for all u e R. If A and B are not 
both 0, we put X(u) — X(u) + B ? iA and direct calculations yield that 
the curve X satisfies equation (j4]). But, as we saw before, that means 
X screw-dilates around and hence X screw-dilates around the point 
— B ? iA ■ So when we have rotation or dilation (or both), introducing 
the translation term doesn't give us any new solutions. 

However, if both A and B are we can take / = 0, g = 1, H = Ct 
and then it is clear that equation (T5]) is satisfied for all (u,t) e R x I. 
Thus, the curve X moves by translation in the direction C under the 
CSF. Hence, the curve is either a straight line with C as its direction 
vector (making the translation vacuous) or a part of the curve can be 
written as a graph of a function defined on a line perpendicular to the 
direction C. However, in [E], it is shown that in this case the curve 
is the (possibly rescaled) graph of the function x i— > — log(cosx), the 
so-called Grim Reaper curve; see Figure [TJ 

So the problem has been reduced to finding all immersed curves 
X : R — ?• C satisfying equation fll]). From now on we will parametrize 
the curves by arc length. For every curve X, we have 

j-(X,T) = l + k(X,N) 

±(X,N) = -h(X,T). 

as 

Thus, if we define the functions 

x = A(X } T) + B(X,N) 
y = -B(X,T) + A(X,N), 
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which can also be written as 

x + iy = (A - iB)((X, T) + i(X, N)), 

then they satisfy 

x = ky + A 
y' = —kx — B. 

Note that we always have 

X = {X,T)T+(X,N)N 

= T((X,T)+i(X,N)) 

_ je x + W 
A-iB 1 

where 0{s) = f° k(s')ds' + 9 and T(0) = e ie °. 

Now, we are looking for curves X which satisfy the equation x = k 

and by what was shown above, the following method yields all possible 

curves with that property. We let x and y be the unique solution to 

the system of differential equations 

x = xy + A 

i 2 D 

y = —x — B 

for any initial values x(0) = xq and y(0) = yo- Note that 

d rr, n Ax — By /—- — 

— Jx 2 + y 2 = = < VA 2 + B 2 

so x and y can't blow up in finite time and hence the solution is defined 

on all of R. Then we put 8(s) = J x(s')ds' + 9q for any 6$. Finally, 

we define the curve as 

46 x + l V 



Note that 



X' 



X = e v 

A-iB 

ie iO'(x + iy) + x' + iy' 
6 A~^iB 

ie ix(x + iy) + xy + A + i(—x 2 — B) 

6 A^W 

je 
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so X is parametrized by arc length with tangent T = e %e and hence the 
curvature k is equal to 9' = x. Finally, 

A(X, T) + B(X, N) = (X, (A + iB)T) 

= Re{X{A - iB)e~ ie ) 
= x 
= k 

so equation (j4j) is indeed satisfied. Hence, we have proved the following. 

Theorem 2.1. For each value of A and B there exists an immersed 
curve X satisfying equation (J2J) . 

In fact, unless A and B are both 0, we have a whole family of im- 
mersed curves satisfying equation (J3J). The rest of this paper will be 
devoted to classifying these curves. 

There is a lot of symmetry in the problem. Each curve is determined 
by the three parameters xq, yo and 9q. However, changing 6q only 
rotates the curve so in order to avoid getting identical curves we will 
from now assume 9q = 0. Also note that rescaling the value of (A, B) 
by c corresponds to rescaling the curve X by c~2. 

Define the functions F(x, y) = xy + A and G(x, y) = —x 2 — B such 
that we have 

x = F(x,y) 
y = G(x,y). 

Note that F(—x,—y) = F(x,y) and G(—x,—y) = G(x,y) and there- 
fore s i— > —(x(—s),y(—s)) is also a solution to the system of ODEs, 
which of course just corresponds to parametrizing the curve X back- 
wards. This symmetry will be used to simplify some arguments. 

It will turn out to be convenient to work also with the functions 
{X, T) and (X, N) which from here on will be denoted by r and v 
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respectively. Recall that 

x + iy = (A — iB){r + iv) 
and since the curvature equals x, their derivatives are given by 

t' — 1 + xv 



v = —XT. 



Hence, they are solutions to the system of differential equations 

r' = 1 + Atv + Bv 2 
v' = -At 2 - Btu 

and their initial values satisfy the equation xo + iyo = (A — iB)(ro + ivo). 
The curve X can therefore be written as 

(5) X = e ie ^V- = e w (T + iv) = e l<t> r 

y ' A-iB K J 



where 



ixi = V- 



r z + u z 



\J x 2 + y 2 



y/A 2 + B 2 

and 

<fi = + arg(r + iv) = 9 + arg(x + iy) + arg(A + iB). 



All possible values of A and B will be divided into 4 cases which 
yield curves with significantly different properties and behaviour under 
the flow. Those are 



A ^ and B > 
A ^ and B < 
A = and B < 
A = and B > 



Curves which rotate and expand 
Curves which rotate and shrink 
Shrinking curves 
Expanding curves 



The case A = B = only yields straight lines and is therefore omitted. 

Before we look at each case separately, let's address a simple ques- 
tion: Can the curves have double points? Note that a double point 
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would have to be transversal by the uniqueness of our system of ODEs. 
Now, assume X is any curve with a transversal double point. Then 
X has a loop, I\ Let the curve be parametrized such that the region 
enclosed by the loop is on the left side of the curve and let a G (0, 2ti) 
be the interior angle of the loop at the double point. If we let the curve 
X flow under the CSF then the time derivative at time of the area 
enclosed by the loop is 



(6) - kds = -(tt + o) < 0. 

Hence, the area enclosed by any loop is always strictly decreasing under 
the CSF. That means that if a curve expands under the flow (corre- 
sponding to B > 0), then it can't have any double points since the area 
enclosed by the corresponding loops would be increasing. 



3. Curves which rotate and expand 

For each A^O and B > 0, we have the following description of the 
curves satisfying equation (TJJ. 

Theorem 3.1. The curves are embedded, have one point closest to 
the origin and consist of two arms coming out from this point which 
strictly go away from the origin to infinity. Each arm has infinite total 
curvature and spirals infinitely many circles around the origin. The 
curvature goes to along each arm and their limiting growing direction 
is B + iA times the location. 

The curves form a one- dimensional family parametrized by their dis- 
tance to the origin. 

If B = then under the CSF the curves rotate forever with constant 
angular speed A. 

If B > then under the CSF the curves rotate and expand forever 
with angular function ^ \og(2 Bt + 1) and scaling function \J1Bt + 1. 

Curves of this kind can be seen in Figures IHlfTTl 
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The proof will be the result of a series of lemmas. By reflecting 
the curve if necessary, we can assume A > 0. We start with a simple 
observation. 

Lemma 3.2. The curvature x has at most one zero, is negative before it 
and positive after it. Moreover, x has at most two extrema, a minimum 
below and a maximum above 0. 

Proof. The first statement is true, since if x(s) = then x'(s) = A > 
so x has to grow through each of its zeros and thus can only have one. 
If x'(s) = then x"(s) = —x(s)(x 2 (s) + B) so all extrema above are 
maxima and all extrema below are minima. Therefore, x has at most 
one minimum and at most one maximum which are below and above 
respectively. □ 

Lemma 3.3. The curvature x is bounded. 

Proof. Can x > 0, x' > hold for all s > so for some so? No, since then 
y' < — x 2 (s ) < so y decreases at least linearly and hence x' = xy + A 
will sooner or later reach 0. Similarly, x > 0, x' < cannot hold for 
all s < Sq since then y' < —x 2 (s ) < so y increases at least linearly 
in the backwards direction and hence x' = xy + A will reach for 
a small enough s. From this and Lemma 13.21 we conclude that if x 
reaches a positive value then it has to have a global maximum. By 
symmetry it is clear that if x reaches a negative value then it has a 
global minimum. □ 

Lemma 3.4. \im s ^± OQ y = =poo. 

Proof. If B > 0, then y' < — B < so the statement is obvious. 
Therefore we look at the case B = 0. Since y' = —x 2 < 0, y is 
decreasing. Assume y is bounded from below so that lim^oo y is finite. 
Now, y" = —1x 2 y — 2Ax is bounded, since both x and y are bounded. 
Hence, by Barbalat's lemma (see [SL]) lim^oo y' = 0. But y' = —x 2 
so lim^oo x = 0. However, x' = xy + A so lim^oo x' = A, which 
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contradicts lim^oo x = 0. Hence Hindoo y = — oo and by symmetry, 
lim^.oo y = oo. □ 

Lemma 3.5. lim^-too r = oo and r has exactly one extremum, a global 
minimum. 



Proof. The first statement follows from Lemma [331 since r (A +B ) = 
x 2 + y 2 . To see the second one, note that ^r 2 = 2r and if r(s) = 
then t'(s) = Bu 2 {s) + 1 > 0. □ 

This lemma implies that the curve X has exactly one point closest 
to the origin. The curve consits of two arms coming out from this point 
and each arm is strictly going away from the origin to infinity. Hence, 
each of these arms is embedded. Since we already know that X doesn't 
have double points (because B > 0) we have the following interesting 
result. 

Corollary 3.6. The curve X is embedded. 

Our next goal is to show that the curvature x goes to along each 
arm. 

Lemma 3.7. The curvature x has a negative global minimum, a posi- 
tive global maximum and no other extrema. 

Proof. If x < for all s then x' = xy + A>A>0 for all s large 
enough (by Lemma I3.4[) and hence x will eventually attain the value 
0, a contradiction. Similarly, x > can't hold for all s. Therefore, 
x takes both negative and positive values and the argument from the 
proof of Lemma 13.31 shows that it has a global minimum and a global 
maximum. By Lemma 13.21 these are the only extrema. □ 

Lemma 3.8. lim s _ >= |- 00 v = =poo. 

Proof. Follows from Lemmas 13.31 and 13.41 since v = ^\^. ■ □ 

Lemma 3.9. lim^-too x = 0. 
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Proof. By Lemma 13. 7\ we know that x has a finite limit in each direc- 
tion. If the limit of x when s — > oo is not 0, then x > 5 > for all 
s > Sq. By Lemma [3. 8 \ we can pick so such that v < — | for all s > So- 
But then r' = 1 + xz/ < 1 — x| < — 1 for s > s so lim s _ >00 r = — oo. 
But x = At + Bz/ so this yields lim^oo x = — oo, a contradiction. The 
other case follows by symmetry. □ 

Now we can determine the behaviour of the two arms far away from 
the origin. First let us state a basic limit result. 

Lemma 3.10. lini^-too ^ — 

um s _> ±oc - - it VA2+B -, 

]\ rn v_ L —A 



A 



V 



Proof. Follows from Lemmas 13.41 and 13.91 since r = A f +B f and 

Bx+Ay n 

By equation (jSJ) and the previous lemma, we immediately get the 
following corollary which shows the limiting growing direction of the 
arms. 

Corollary 3.11. lim^ ±00 f = ±-/jj^s. 

Lemma 3.12. lims^-too cj) = +oo. 

Proof. Easy calculations give dlo t r) = — ^ which goes to ^ (or +oo 
when B = 0) as s — > ±oo. Hence the result follows from Lemma 
1331 □ 

The above results show that each arm spirals infinitely many circles 
out from the origin and that its limiting growing direction is B + %A 
times the location. Since <fi = 9 + arg(r + iv) and we just saw that 
arg(r + iv) has finite limits as s — > ±oo we get the following corollary. 

Corollary 3.13. lim s ^± OQ 9 = +oo. 

In other words f°°kds = +oo and f s ° k ds = — oo so each arm has 

J so J— oo 

infinite total curvature. 
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This result concludes the proof of the theorem. Note that in the 
purely rotational case (B = 0) simple calculations yield that the func- 
tion r 2 — \Q is constant, which is an interesting fact. This also gives 
another proof of the fact that the curve can't cross itself. 



4. Curves which rotate and shrink 

For each A ^ and B < 0, we have the following description of the 
curves satisfying equation 01]). We omit the the circle of radius -7=??, 
where the rotation is vacuous. 

Theorem 4.1. In this case there are two types of curves. 

1) Curves such that the limiting behaviour when going along the curve 
in each direction is wrapping around the circle of radius -/=m, clockwise 
if A < and counterclockwise if A > 0. These curves form a one- 
dimensional family. 

2) Curves such that the curvature never changes sign and the two 
ends behave very differently. One end wraps around the circle with 
radius -7= in its limiting behaviour, clockwise if A < and counter- 
clockwise if A > 0. The other end spirals infinitely many circles around 
the origin out to infinity and has infinite total curvature. The curva- 
ture goes to along it and its limiting growing direction is —B — iA 
times the location. There is at least one curve of this type and we call 
it the comet spiral. 

Under the CSF these curves rotate and shrink with angular function 



_A_ 
2B 



log(2£>£ + l) and scaling function \/2Bt + 1. A singularity forms at 
time t = — ^g. The curves of type 1 are bounded so they disappear into 
the origin. 

Curves of this kind can be seen in Figures [1211191 



The proof will be given through a series of lemmas. By reflecting 
the curve X if necessary, we may assume A > 0. By rescaling, we may 
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assume A = 1 and B = —(3 2 for some (3 > 0. Then we have the system 

x' = xy + 1 

y = —x + p . 

By symmetry, it will suffice to look at the behaviour of the solutions in 
the right half-plane. There the system has a fixed point (/?, — 4) which 
corresponds to X being the circle of radius 4 parametrized counter- 
clockwise (making the rotation vacuous). The eigenvalues of the lin- 
earization of the system around the fixed point are 



-1 ± y/1 - 8/3 4 
W ' 

Since the real part is always negative, the fixed point is a sink, i.e. asymp- 
totically stable. Moreover, it is a spiral if [3 > 8" 1 / 4 and a node if 
< (3 < 8 -1 / 4 . It will turn out to be asymptotically stable in the 
whole right half-plane. We begin with a simple lemma. 

Lemma 4.2. The curvature x has at most one zero, is negative before 
it and positive after it. Moreover, if x'(s) = then x has a maximum 
at s if x(s) > (3 and a minimum if < x(s) < (3. 

Proof. The first statement is clear since x'{s) = 1 if x(s) = 0. The 
second statement follows from the fact that x"(s) = —x(s)(x 2 (s) — (3 2 ) 
if x'(s) = 0. □ 

Lemma 4.3. If xo = 0, then lim s _ >00 (x,y) = (p\— i). 

Proof. First assume x < j3 always holds. Then y' > so y is increasing. 
If y attains the value 0, then x' > 1 which implies x will eventually 
reach the value (3, a contradiction. Hence, y is bounded so it has a finite 
limit. It is clear that y" is bounded since both x and y are bounded 
and hence, by Barbalat's lemma y' — > 0. But then x — > f3 and since 
both x and y have limits, the limit has to be the fixed point (/3, — 4). 
In the other case, x reaches (3. Let s\ be the first time that happens. 
Then 2/'( s i) = anc i we m ust have y{s-y) > —4 since x'(si) > 0. Again 
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we have two possibilities. First assume x > (5 holds for all s > s±. 
Then y' < so y is decreasing. If y reaches —4 then we will have 
x' < —5 < so x goes down to (3, a contradiction. Therefore y is 
bounded from below and has a finite limit. Also note that x can't 
always be increasing since then y' would be decreasing and y would 
go to — oo. Therefore, x reaches a maximum and after that will be 
decreasing, by Lemma 14.21 Hence, x also has a finite limit and thus 
(x,y) has the fixed point as its limit. 

The other possibility is that x reaches (3 again. Let S2 be the first 
time that happens. Then y'(s2) = and we must have y(s2) < — \ 
since x\s-i) < 0. By repeating the argument from the first paragraph 
we get that either (x, y) goes directly to the fixed point or there exists 
a S3 such that x(ss) = (3 and y(ss) > —-s. Also note that since x = 0, 
x > for all s > and the (x, y) trajectory can't cross itself, we must 
have 2/(53) < y(si), i.e. the curve spirals inwards. 

Now, let K be the compact region in the plane with boundary con- 
sisting of the (x, y) trajectory between s\ and S3 and the line segment 
between those two endpoints. Then it is clear that the trajectory after 
S3 will never leave K since it can neither cross its own path nor the 
line segment because there the arrows of the (F, G) vector field are 
all pointing into K. Hence the result of the lemma follows from the 
Poincare-Bendixson theorem (see [GHj ) and the following lemma. □ 

Lemma 4.4. There are no periodic orbits in the right half-plane. 

Proof. Define the scalar field H(x, y) — - which is smooth in the right 
half-plane. Since 



d d 1 

dx dy x 2 



does not change sign in the right half-plane, by Dulac's criterion (see 
[GHJ) there can be no periodic orbits there. □ 
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Fix the trajectory with xq = yo = 0, the symmetric trajectory. How 
do the other trajectories in the right half-plane behave? The following 
lemma gives the answer when we go forward in time. 

Lemma 4.5. The fixed point is asymptotically stable in the whole right 
half-plane. In other words, if xq > then lim s _ s>00 (x, y) = (/?, —4). 

Proof. In Lemma [4731 we proved this result for trajectories in the right 
half-plane which cross the y-axis. However, the only place we used that 
extra condition was to show that the trajectory can't spiral outwards. 
But if one trajectory spirals inwards, then all the others also have to 
spiral inwards since otherwise they would be forced to cross the one 
spiraling inwards. Hence the same argument as before works for all 
trajectories in the right half-plane. □ 

But what do the other trajectories do when we go backwards in time? 
Let Q be the region in the plane defined by Q = {(x,y) : < x < /3}. 
By the following lemma, it is clear that if a trajectory leaves Q through 
the right boundary, then it has to come back in. 

Lemma 4.6. If x(s ) > (3, there exists s < s such that x(s) = (3. 

Proof. Assume not. Then x > (3 for all s < Sq so y' < and thus 
y is decreasing (increasing in backwards direction). If y manages to 
reach when going backwards then x' > 1 so x will also reach /3, a 
contradiction. Therefore, y is bounded from above and has a finite 
limit at — oo. If x' < for all s < sq then x is decreasing so y' is 
decreasing and hence y goes to +oo when s —¥ — oo, a contradiction. 
Hence, x reaches a maximum and before that is increasing, by Lemma 
14.21 Thus x also has a finite limit when s — > — oo. However, the limit 
of (x,y) when s — > — oo would have to be a fixed point and the only 
one in the right half-plane is (/?, — i). This can never happen since the 
fixed point repels trajectories when we go backwards in time. □ 

It is also clear that each trajectory will ultimately stop leaving Q 
through the right boundary, since otherwise it would have to cross the 
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symmetric trajectory. After that there are two possible behaviours: 

1) The trajectory leaves Q through the left boundary, i.e. crosses the y- 
axis. By Lemma H~3l and symmetry, it is then clear that lim s _ > _ 00 (a;, y) = 
(— (3, 4), i.e. the trajectory goes to the fixed point in the left half-plane 
which is a source. Since r 2 (l + /3 4 ) = (x 2 + y 2 ), we have 

lim r = — . 

s-»±oo p 

Hence, each arm of the curve X has as its limiting behaviour wrapping 
counterclockwise around the circle of radius 4. These curves form a 
one- dimensional family that can for example be parametrized by the 
value yo when x = 0. 

2) The trajectory stays in Q forever. Then by the following lemma 
we have 

lim x = 0, lim y = — oo and hence lim r = oo. 

s— >— oo s— >— oo s— >— oo 

Hence, the curve X has one arm wrapping counterclockwise around the 
circle of radius 4 and the other one going out to infinity. The same 
arguments as in Lemmas I3.10[ 13.121 and Corollaries 13.111 and 13.131 show 
that the arm has infinite total curvature, spirals infinitely many circles 
around the origin and its limiting growing direction is —B — %A times 
the location. The curvature of X never changes sign. 

Lemma 4.7. The trajectories of type 2 satisfy 

lim x = and lim y = — oo. 

s— >— oo s— >— oo 

Proof. By Lemma 14.21 x has at most one extrema (after entering Q for 
the last time) and since it is bounded, it has a finite limit as s — > — oo. 
Since the curve remains in the area Q, we have y' > so when going 
backwards in time y is decreasing. If y were bounded from below it 
would have a finite limit as s — > — oo. But then (x, y) would have to 
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go to the fixed point, a contradiction since it repels trajectories when 

we go backwards in time. Hence, it is clear that y — > — oo as s — > — oo. 

Now, since x' = xy + 1 it is clear that if linis^.oo x > then 

lim^.oo x' = — oo, a contraditction. Hence, lim^^.oo x — 0. □ 

So the remaining question is: Does there exist a curve of type 2? 

Lemma 4.8. There exists a curve of type 2. 

Proof. The trajectories of type 1 can be divided into two sets, those 
with y > and those with y < (where we assume xq = 0). Since 
all the trajectories in the right half-plane go to the fixed point, the two 
sets will somehow have to meet. The symmetric curve (y^ = 0) gives 
one of their boundaries but there must also be a curve (at least one) 
on their other boundary, i.e. when y goes to oo and — oo. However, 
this curve can be in neither of the two sets since they both have the 
structure of an open interval (||/o| can always be made bigger). Hence, 
the curve does not cross the y-axis and therefore must be of type 2. □ 

This finishes the proof of the theorem. Note that these curves are 
usually not embedded and have loops. By equation ([6]) and the fact 
that the motion is self-similar we get that the area inside each loop 
decreases with the fixed rate — (tc + a) where a G (0, 2-k) is the interior 
angle at the double point. The loop shrinks to a point at time t = — -^ 
and therefore we have the following relationship between the interior 
angle and the area inside the loop 

Area inside loop = . 

2B 

Of course, we only get a = n (smooth loop) when the curve X is the 
circle of radius / . Since a > it follows that each loop has at least 
half the area of the circle, which is an interesting fact. 



18 hoeskuldur p. halldorsson 

5. Shrinking curves 

The remaining two cases are already more or less known (see e.g. |ALj . 
[I] and [U]) but we cover them here for the sake of completeness. 

For A = and B < 0, we have the following description of the curves 
satisfying equation fll]). We omit the straight lines through the origin 
where the scaling is vacuous. 

Theorem 5.1. Each of the curves is contained in an annulus around 
the origin and consists of a series of identical excursions between the 
two boundaries of the annulus. The curvature is an increasing function 
of the radius and never changes sign. 

The curves form a one- dimensional family parametrized by their dis- 
tance to the origin (inner radius of annulus) and are divided into two 
sets: 

1) Closed curves, i.e. immersed^ 1 (Abresch-Langer curves). In addi- 
tion to the circle, there is a curve with rotation number p which touches 
each boundary of the annulus q times for each pair of mutually prime 
positive integers p, q such that | < - < &. 

2) Curves whose image is dense in the annulus. 

Under the CSF these curves shrink with scaling function g(t) = 



\j2Bt + 1 until they disappear into the origin at time t 
Curves of this kind can be seen in Figures HUH 



i 

2B- 



The proof is contained in a series of lemmas. By rescaling if neces- 
sary, we may assume B — — 1. Then we have the following system for 
x and y 

x' = xy 

y' = -x 2 + 1. 

From the first equation we immediately conclude 

(7) x{s) = x e^ y{s ' )ds ' 
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so the curvature x does not change sign. When xq = we just get the 
solution (x, y) = (0, yo+s) which corresponds to X being a straight line, 
making the scaling vacuous. By parametrizing the curves backwards 
if necessary we may assume xq > and hence x > 0. Note that the 
system has the fixed point (1,0) corresponding to X being the unit 
circle parametrized counterclockwise. 
We start with a technical lemma. 

Lemma 5.2. Ifyo > 0, then there exist s\ < < S2 such that y(si) = 
and y(s 2 ) = 0. 

Proof. Assume y > for all s > 0. Then x' > so x is increasing. If 
x < 1 holds for all s > then y' = — x 2 + 1 > so y is increasing. 
But then by equation (J7J), we have x — > oo when s — > oo, contradicting 
x < 1. Hence, x must reach values above 1 which implies y' < —5 < 
so y has to go down to 0, a contradiction. The other case is treated 
similarly. □ 

Lemma 5.3. The solution (x,y) is periodic. 

Proof. Since F(x,—y) = —F(x,y) and G(x,—y) = G(x,y) it is clear 
that s I—)- (x(—s),—y(—s)) is also a solution to the system of ODEs. 
Hence, if y = 0, the initial values are the same so by uniqueness, 
we get that (x(—s),—y(—s)) = (x(s),y(s)). But this means that the 
trajectory is symmetric with respect to reflection across the x-axis. By 
Lemma 15.21 each trajectory crosses the x-axis twice and therefore has 
to be closed, i.e. periodic. □ 

Corollary 5.4. Both the curvature x and the radius r are periodic. 

The phase portrait of our dynamical system therefore consists of the 
y-axis, the fixed point (1, 0) and closed trajectories around (1,0). Note 
that for each of these the curvature x takes its maximum and minimum 
on the x-axis, since x' = xy. 

Direct calculations give us the following lemma. 
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2 . 2 

g +1/ 



Lemma 5.5. The function x e 2 is constant. 

Remember that (A, 5) = (0,-1) so we have r 2 = x 2 + y 2 . Now, 
imagine we parametrize our solution such that yo = 0. By Lemma [5.51 
we have the following relationship between the curvature x and the 
radius r 



2 2 



x = x e 2 

where Xq can be either the minimum of the curvature, x m i n , or its 
maximum, x max . In particular, r is an increasing function of x so we 
have the following corollary. 

Corollary 5.6. The radius takes its maximum and minimum at the 
same time as the curvature x. Moreover, r min = x min and r max = x max . 

So the curve X consists of infinitely many identical excursions be- 
tween the two boundaries of the annulus with inner radius r min and 
outer radius r max . The curvature is an increasing function of the ra- 
dius and reaches its lowest and highest value at the innermost and 
outermost points respectively. 

Now let's examine a whole excursion, i.e. we go from the outer bound- 
ary, touch the inner bounday and go back out the outer boundary. This 
corresponds to the (x, y) trajectory going one period. Therefore, the 
difference between the values of the location angle (p is the same as the 
difference between the values of the tangent angle 9. This angle differ- 
ence, A9, is given by the integral of the curvature x over one period of 
the (x, y) trajectory. Now we have two possibilities: 

1) The angle difference A9 is of the form — for some positive in- 
tegers p, q. Then after q excursions the curve X closes up. Thus, X is 
a closed curve, an immersed S 1 . The number p is the rotation number 
of the curve and q is the number of times it touches each boundary of 
the annulus. 
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2) The angle difference A8 is not of this form. In this case the curve 
never closes up and takes infinitely many excursions. It is dense in the 
annulus. 

In |ALj . Abresch and Langer investigated A9 as a function of r max 
and showed that it is a continuous decreasing function with limits 

lim A6 = y/2n and lim A6 = n. 

Therefore, we get a closed curve for all mutually prime positive integers 
p, q such that | < - < *2.. This finishes the proof of the theorem. 



6. Expanding curves 

For A = and B > 0, we have the following description of the curves 
satisfying equation (jlj). 

Theorem 6.1. Each of the curves is convex, embedded and asymptotic 
to the boundary of a cone with vertex at the origin. It is the graph of 
an even function. 

The curves form a one-dimensional family parametrized by their dis- 
tance to the origin. 

Under the CSF these curves expand forever as governed by the scaling 



function g(t) = \?2Bt + 1. 

Curves of this kind can be seen in Figures [2] and [3j 

As before, the proof will be given through a series of lemmas. By 
rescaling we may assume B = 1. Then (x, y) is a solution to the system 

x = xy 

I 2 i 

y = —x — 1. 
From the first equation we immediately get 

(8) x{s) = x e^ y{s ' )ds ' 
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so the curvature x does not change sign. When xq = we just get the 
solution (x, y) = (0, yo — s) which corresponds to X being a straight line 
(making the scaling outwards vacuous). By parametrizing the curve X 
backwards if necessary, we may assume x > and hence x > 0. 

Since y' < — 1, it is clear that \im s ^± OQ y = =)=oo so by (jHl) we get 
lim^-i-oo x = 0. Also, since -F(x, — y) = —F(x, y) and G(x, — j/) = 
G(x,y), we get the following as in the proof of Lemma [5731 if we pa- 
rametrize the solution such that y = 0. 

Lemma 6.2. The curvature x is an even function and y is an odd 
function. 

Thus, the angle 9 is an odd function (remember we assume 9q = 0) 
and since X is given by (|5l), we have the following corollary. 

Corollary 6.3. The curve X is symmetric with respect to reflection 
across the y-axis. 

Direct calculations yield this lemma. 

x -\-y 

Lemma 6.4. The function x e a is constant. 

Since (A, B) = (0,1) we have r 2 = x 2 + y 2 . Thus, the following 
relationship between the curvature x and the radius r holds 



x = x e 2 

where Xq = x max , the maximum value of the curvature, since x' = xy. 
In particular, r is a decreasing function of x so the following is true. 

Corollary 6.5. The radius r takes its minimum when the curvature x 
takes its maximum. Moreover, r min = x max . 

Lemma 6.6. The total curvature of X is strictly less than n. 
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Proof. The total curvature is given by the integral 



x(s)ds = 2 
<2 
= 2 
= 2 

= 7T. 
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where we put u = ye 1 



min 



1 in the second to last equality. □ 

This result and Corollary 16.31 yield that X is the graph of an even 
function and hence embedded. It also follows that X is convex. In [I], 
Ishimura shows that the curve is asymptotic to a cone with vertex at 
the origin, finishing the proof of our theorem. Moreover, he shows that 
the total curvature is an increasing continuous function of the distance 
to the origin which gives a homeomorphism between [0, 00) and [0, 71"). 



7. Figures 




-10 12 3 



Figure 1. The Grim Reaper curve, which translates under 
the flow. 
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Figure 2. A = 0, B = 1 (expands) The curve is asymptotic 
to the boundary of a cone with vertex at the origin. 
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Figure 3. A = 0, B = 1 (expands) The total curvature 
increases with the distance to the origin. 
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Figure 4. A = 0, B = —1 (shrinks) Abresch-Langer curve 
with p = 2 and q = 3. 
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Figure 5. A = 0, B = —1 (shrinks) Abresch-Langer curve 
with p = 7 and q = 10. 
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Figure 6. A = 0, B = —1 (shrinks) Abresch-Langer curve 
with p = 20 and q = 31. 
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Figure 7. A = 0, B = — 1 (shrinks) Abresch-Langer curve 
where the values of p and q are very high. 



SELF-SIMILAR SOLUTIONS TO THE CURVE SHORTENING FLOW 27 




Figure 8. A = 1, B = (rotates) The symmetric curve is 
Altschuler's yin-yang spiral. 




Figure 9. A = 1, B = (rotates) As the distance to the 
origin increases, the tip of the spiral looks more and more 
like the Grim Reaper. 
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Figure 10. A = 1, B = 0.25 (rotates and expands) The 
symmetric curve looks like a less dense version of the the 
yin-yang spiral. 
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Figure 11. A = 1, B = 0.25 (rotates and expands) It gets 
thinner as the distance to the origin increases. 
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Figure 12. A = 1, B = -0.05 (rotates and shrinks) The 
symmetric curve is similar to the yin-yang spiral but only 
spirals out to the limiting circle. It looks embedded for values 
of B close to 0. 




Figure 13. A = 1, B = —0.05 (rotates and shrinks) Here we 
start on the outside of the circle and spiral inwards towards 
it. 
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Figure 14. A = 1, B = —0.05 (rotates and shrinks) This 
is the comet spiral. It seems to be embedded for values of B 
close to 0. 
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Figure 15. A = 1, B = —0.25 (rotates and shrinks) Here 
the symmetric curve is similar to the yin-yang symbol. Is it 
still embedded? 
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Figure 16. A = 1, B = —0.25 (rotates and shrinks) A 
typical curve of this kind looks like a spirit coming out of a 
circle. Each loop has at least half the area of the circle. 
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Figure 17. A = 1, B = — 0.25 (rotates and shrinks) Is the 
comet spiral still embedded here? 
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Figure 18. A = 1, B = — 5 (rotates and shrinks) The sym- 
metric curve is no longer embedded when B is small. 




1.5 



Figure 19. A = 1, B = — 5 (rotates and shrinks) Similarly, 
the comet spiral is not embedded when B is small. 
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